UT-11-38 
TU-894 
IPMU-11-0188 



Exploring Supersymmetric Model 

O ■ with Very Light Gravitino at the LHC 

(N ■ 

> 
O 



Masaki Asano^'''^), Takumi Ito^''-'^), Shigeki Matsumoto^'^) 
and Takeo Moroi^^ '^) 



^"•^11. Institute for Theoretical Physics, University of Hamburg, 
Luruper Chausse 149, DE-22761 Hamburg, Germany 
^ , ^^"^ Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 

^ ■ Department of Physics, Tohoku University, Sendai 980-8578, Japan 

; ^'^HPMU, TODIAS, University of Tokyo, Kashiwa, 277-8583, Japan 



^ ' Abstract 



The low-scale gauge mediation scenario of supersymmetry breaking pre- 
dicts very light gravitino, which makes the next lightest supersymmetric par- 
ticle (NLSP) quasi stable. We study the LHC phenomenology of the case 
that the NLSP is the stau. When the mass of the gravitino is of the order of 
lOeV, the decay length of stau is about 0.1-1 mm, so that it decays before 
reaching the first layer of the inner silicon detector. We show, however, that, 
with utilizing the impact parameter of r-jets from the stau decay, it is possi- 
ble to determine the mass spectrum of sparticles precisely. It is also possible 
to estimate the lifetime of the stau by observing distribution of the impact 
parameter. 



1 Introduction 



Large Hadron Collider experiment (LHC) is now operating and reports many im- 
portant results on new physics beyond the standard model (SM). Although positive 
signals have not been reported so far, those are expected to be found in near fu- 
ture, because the hierarchy problem of the SM strongly suggests the existence of 
new physics at the TeV scale or below. On the other hand, many new physics mod- 
els have been theoretically proposed. Among those, the supersymmetric model is 
very attractive because it guarantees the stability of the Higgs mass to its radiative 
corrections and gives a clue to solve the hierarchy problem. In addition, the su- 
persymmetry (SUSY) plays a crucial role to realize the grand unification of known 
gauge interactions of the SM at a certain high energy scale. 

Details of supersymmetric model, such as the mass spectrum of sparticles, depend 
highly on how SUSY is broken. So far, a variety of SUSY breaking mechanisms 
has been proposed [I]. Among those, the gauge mediation scenario [2] attracts 
an attention, because it gives a solution to dangerous SUSY flavor problems. In 
this scenario, the breaking occurs at lower energy scale than those of other SUSY 
breaking scenarios, so that the superpartner of graviton, the gravitino, is likely to be 
the lightest supersjTiimetric particle (LSP). The gravitino mass is predicted to be in 
the range between O{10)eV and 0{l)GeV. In this article, we focus on the low-scale 
gauge mediation model providing a gravitino with O{10)eV mass. Such a scenario is 
well motivated because it is completely free from severe cosmological constraints [3] 
such as Big-Bang Nucleosynthesis [1] and large scale structure formation of our 
universe [5]. 

Collider signals of the low-scale gauge mediation scenario depend on what the 
next lightest superparticle (NLSP) is, which decays only into gravitino and its super- 
partner. Though there are many candidates for NLSP, we focus on the stau NLSP in 
this article, which is predicted in wide parameter region of the scenario. When the 
gravitino mass is of (9(10)eV, the stau NLSP decays into a r-lepton and a gravitino 
with the lifetime of 10^^^-lO^^^sec. The decay length (the lifetime times the speed 
of light) of the stau NLSP is therefore much shorter than the typical size of collider 
detectors, and the traditional supersymmetric signal, namely, multi-jets associated 
with missing energy and r-leptons, is expected at the LHC experiment. Such a signal 
is, however, generally predicted in various SUSY breaking scenarios. 



1 



We show in this article that, even if the NLSP decays before reaching inner track- 
ers of colhder detectors, we can use the impact-parameter information about the 
decay products of the NLSP to study various properties of superparticles. In partic- 
ular, the impact parameter is available for charged tracks caused by decay products 
of r-lepton at the stau NLSP decay. If the decay product of stau is found to have 
large impact parameter, it strongly suggests that the underlying SUSY breaking 
scenario is low-energy gauge mediation. Furthermore, the impact parameter is also 
utilized to precisely measure the spectrum of sparticles such as squark, neutralino, 
and stau masses. This is because two tau leptons produced by the cascade decay of a 
squark can be distinguished with each other by using the impact parameter. In addi- 
tion, we may be able to determine the lifetime of the NLSP (i.e., stau in the present 
study) using the impact parameter distribution. When mass and lifetime of the stau 
NLSP are measured, it is possible to determine the gravitino mass assuming that the 
stau decays into gravitino and tau. The scale of SUSY breaking in the low-energy 
gauge mediation scenario is, therefore, obtained. It has been already shown that 
such studies can be easily performed once the e+e~ linear collider becomes available 
[6]. Here, we consider the case of the LHC We will see that the measurement of 
the mass spectrum as well as the determination of the lifetime of the NLSP can be 
performed at the LHC with the help of impact parameter information. 

This article is organized as follows. In the next section, we consider some prop- 
erties of the stau NLSP in the low-energy gauge mediation scenario and discuss how 
the impact parameter from the NLSP decay is utilized in determinations of sparticle 
masses and NLSP lifetime. Our simulation framework is summarized in section [3l in 
which a representative point and several strategies to reduce backgrounds are shown. 
In section HJ simulation results for the measurements of sparticle masses and lifetime 
of the stau NLSP are discussed. Section |5] is devoted to summary of our studies. 

2 Utilizing impact parameter 

In this section, we discuss how the impact parameter is utilized in order to determine 
the mass spectrum of sparticles and the lifetime of stau NLSP. We first briefly review 
some properties of the stau NLSP and define the impact parameter. Then, we discuss 
basic strategies for the measurement of the mass spectrum and the lifetime of NLSP 
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with the use of the impact parameter. 



2.1 Stau NLSP and impact parameter 

Since the LHC is a hadron colhder, colored sparticles such as squarks and gluino are 
expected to be produced at first, which decay into the stau NLSP, the super partner 
of r-lepton, through several cascade channels. The stau NLSP then decays into a 
r-lepton and a gravitino with the following lifetime, 



where Mpi ~ 2.4 x 10 GeV, , and m^i2 are reduced Planck mass, stau mass, and 
gravitino mass, respectively. It turns out from above formula that the decay length 
(the lifetime x the speed of light) of the stau NLSP is estimated to be ~ (9(100)/im 
when gravitino and stau masses are ~ lOeV and ~ lOOGeV, respectively. On the 
other hand, the decay length of r-lepton which is one of main backgrounds against 
the stau signal, is 87/im, so that the decay of stau NLSP into very light gravitino 
can be, in principal, detected if we can reduce SM backgrounds efficiently. 

Since the decay length of the stau NLSP is, at most, (9(l)mm in the parameter 
region of our interest, the stau NLSP decays before reaching the first pixel detector, 
which is located at 5cm (4cm) away from the beam line in the ATLAS detector [7| 
(CMS detector [8J). The lifetime of the stau NLSP is, as a result, difficult to be mea- 
sured using methods usually applied to detect long-lived particles, such as methods 
by observing charged tracks [HI [IDl [HI 112] ■ On the other hand, the lifetime of the 
stau NLSP may still be possible to be determined using the distribution of the im- 
pact parameter, which is obtained by r-jets from the stau NLSP decay. The impact 
parameter is defined as the shortest distance to the track from the interaction point. 
The positional resolution of the ATLAS detector along the longitudinal direction is 
~ 100/im, which is not good compared to that along the transverse direction, 
Ar ~ 10/xml3 Thus, we use the transverse impact parameter which is defined by 



r, = AM Ml, 




(1) 



di = — 




^ Details of those performances are found in the section " Tracking" in Ref . [7] . 
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Figure 1: Typical decay chain in the low-scale gauge mediation scenario. 

where x;^ and are transverse decay point of the /-th r-lepton and transverse 
momentum of the tau-jets from the J-th r-lepton decay, respectively. Note that the 
summation over the index / should not be taken here. The above formula is used in 
our simulation studies, which will be presented in following sections. We expect that, 
at the LHC experiment, the distribution of the impact parameter dj is obtained by 
measuring the shortest distance (projected onto the transverse-plane) to the r-jet 
track from the interaction point. 

2.2 Impact parameter for mass measurements 

We next consider how the impact parameter is utilized in mass measurements of 
sparticles. At the LHC, colored sparticles such as gluino and squarks are expected to 
be produced copiously, and non-colored sparticles are then produced through cascade 
decays of the colored ones. The chain of the cascade decay is, for example, composed 
of following processes; First colored sparticle decays into a neutralino/chargino by 
emitting a quark which is observed as a jet. Next a neutralino/chargino decays into a 
slepton by emitting a lepton. Finally, a slepton decays into a LSP by again emitting 
a lepton. It is needless to say that the LSP passes through the detector without 
giving any signatures, which is, instead, observed as a missing energy. This chain 
(called "golden mode") is frequently used to measure the mass spectrum of sparticles 
in various supersymmetric scenarios by using several kinematical endpoints [T^ [15]. 

In the case of the low-scale gauge mediation scenario, we have a similar decay 
chain. One of the examples is shown in Fig. [H where qr, x\i ^"^^ G are right- 
handed squark, lightest neutralino, lightest stau, and gravitino, respectively. The 
character j denotes a jet which originates in a quark from the qr decayS Following 

^In the gluino production event, the gluino decays into a squark by also emitting a jet. Because 
the mass difference between gluino and squark is much smaller than that between squark and 
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the terminology used in studies of the golden mode, we call r from the x^-decaj 
^(near) ^^^^ that from the fi-decay t^^^^\ Importantly, the impact parameter of t^^^^^ 
is expected to be larger than that of r^"^'^^"^\ which is of great help for the event 
reconstruction. 

In our analysis, we focus on the signal from the decay of a squark shown in 
Fig. [H Kinematics of its decay chain result, governed by following four 

sparticle masses; the masses of squark (m^), lightest neutralino (^^^o), lightest stau 
(m^J, and gravitino {m^/2). Since the gravitino mass is of the order of lOeV, only 
the upper bound on the mass is expected to be obtained. On the other hand, the ex- 
istence of large impact parameters in signal events strongly suggests the (low-scale) 
gauge mediation scenario. We therefore perform our analysis with simply postulat- 
ing that the gravitino mass is much smaller than those of other sparticles, namely, 
with treating the gravitino as a massless particle. Three independent kinematical 
endpoints are then enough to determine the mass spectrum of sparticles. With the 
help of the impact parameter, many kinematical variables are now available. Among 
those, we use the invariant mass between two r-leptons (M^( near) ^(far) ), that between 
jet and near r-lepton (AfjY(near)), and the Mt2 variable from leading two jets {MT2,jj)- 

The upper limit on the invariant mass M^(noar)^(far) is given by 

M^"?ra.),(fao = m^o^l - m?ym|o, (3) 

where the mass of r-lepton is set to be zero in above formula. In addition, the upper 
limit on the distribution of the invariant mass MjY(ncar) is given by the following 
formula, 

M-r^ea.) = m,y(l-m|o/m?) (l-mym|o). (4) 

The last kinematical variable used in our analysis is the Mt2 variable [16] constructed 
from highest two jets; denoting the momenta of highest two jets as p and p', we define 

MT2,jj{m^iss) = min [max {Mr(pT, kr), Mr(pr, k^)}] , (5) 

where Mt is the transverse mass and mmiss is the "test mass." Here, because we 
construct the Mt2 variable only from highest two jets, should be understood as 

neutralino in the parameter region of our interest, we can discriminate between a (soft) jet from 
the ghiino decay and a (hard) jet from the squark decay. 
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the vector sum of transverse momenta of all the activities other than highest two 
jets and missing momentum ^t'- V't = P^t + Zli Pr^jet i + Yl,i P'r-jet %■ The upper 
limit of this variable is then given by 



-^r2,jj(^miss) 



2m^ 



+ 



Ira, 




+ mi 



(6) 



Using three kinematical endpoints given in eqs.(l3]), (jlj) and (E]), we fit the spar- 
ticle masses m^, "n^^o, and vtif^. We will see that, though the cascade chain always 
involves rs as lepton emissions, the spectrum can be determined accurately because 
of information about the impact parameter. 



2.3 Impact parameter for stau lifetime measurement 

Information about the lifetime of the stau NLSP is imprinted in the distribution of 
the impact parameter. The impact parameter, however, depends not only on the 
lifetime of the stau NLSP but also on its mass and velocity. With the use of the 
strategy discussed in previous subsection, the stau mass is measured precisely. On 
the other hand, since the gravitino produced from the stau decay cannot be detected, 
the velocity cannot be determined on event-by-event basis, which makes it difficult 
to determine the lifetime using the impact parameter distribution. 

Once the mass spectrum of the superparticles are known, however, we expect 
to acquire information about the velocity distribution of the f\ in supersymmetric 
events. In the present case (where the mass spectrum of a simple gauge mediation 
model is assumed), we may understand that the underlying scenario is indeed the 
low-scale gauge mediation from the experimentally measured mass spectrum as well 
as the confirmation of the existence of long-lived stau. Even if we cannot specify 
the complete structure of the underlying model, we may still be able to measure 
the masses of superparticles which are most important for the determination of the 
velocity distribution of t\ (i.e., the masses of g, Xi? and fi), as we have discussed in 
the previous section. Then, once those information becomes available, one will be 
able to obtain the velocity distribution with, for example, Monte Carlo simulation. 

In our analysis, we assume that the velocity distribution of f\ can be understood 
once the superparticles are discovered. The detailed study of the methods of de- 
termining the velocity distribution is beyond the scope of this paper, so we simply 
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assume that the averaged velocity of the produced fi can be obtained with some ac- 
curacy and determine the Hfetime using the fi. Although a better determination of 
the lifetime of fi may be possible if we can somehow obtain and use the information 
about the velocity distribution of fi, we can still have a relatively good determina- 
tion of the lifetime using the averaged velocity as we will describe. The procedure 
to measure the lifetime of the stau NLSP is therefore the following. 

(i) We first assume that the averaged velocity of the stau NLSP, denoted as Pfi , 
is somehow understood. Then, we generate fi with the fixed velocity and 
make signal templates of the distribution of the transverse impact parameter 
dr- The template is prepared for wide range of stau lifetime. The production 
angle of the stau NLSP is assumed to be isotropic in generating the events for 
the template. 

(ii) Distribution of the impact parameter expected at the LHC experiment is ob- 
tained by using Monte Carlo simulation. 

(iii) Comparing the templates obtained in the (i) with the actual distribution ob- 
tained in (ii), we study how well we can constrain the lifetime of fi by x^- 
analysis. By varying the value of used in making the templates, we also 
discuss the uncertainty related to the determination of the velocity distribu- 
tion. 

Here, we have a few comments on the above method. First comment is on the 
effect of gluino production. The mass difference between gluino and squark is much 
smaller than that between squark and neutralino in the parameter region of our 
interest. In addition, it will be possible to select signal events with the desirable 
squark decay chain {q ^ Xi ~^ ^ G) by applying appropriate kinematical cuts. 
Thus, the averaged value of the boost factor is expected to depend weakly on the 
gluino mass. We have checked this statement quantitatively by simulating signal 
events with several choices of the gluino mass. Second comment is on how the 
transverse impact parameter dj depends on the production angle of the stau NLSP. 
One might worry if we may compare the actual impact parameter distribution with 
theoretical templates obtained by postulating isotropic distribution of the production 
angle. We have checked that this potential problem can be solved by only using r-jets 
with small pseudo-rapidity. 
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3 Simulation framework 



Before showing our results, we summarize the framework of our simulation study. We 
first mention a representative point and simulation tools used in the study. Next we 
discuss the strategy to suppress combinatorial backgrounds of signal events caused 
by the existence of two decay chains. We finally consider the SM backgrounds and 
discuss kinematical cuts used to reduce those backgrounds. 

3.1 Representative point &; simulation tools 

The representative point used in our simulation study has been chosen by adopting 
the minimal model of the gauge mediation symmetry breaking [2] with the following 
underlying parameters; the SUSY breaking scale {A = 30TeV), the messenger mass 
scale (Mjness = 300TeV), the number of SU(5) messenger fields {N^ = 5), and the 
ratio of vacuum expectation values of two Higgs fields (tan/3 = 15). The ISA JET 
package [17] is used in order to calculate the spectrum and branching fractions of 
sparticles. Resultant masses and branching fractions of sparticles relevant to the 
study are summarized in Table [H for the case of the gravitino mass of 9.7eV which 
corresponds to the decay length of the stau of SOO/xm. It can be seen that the model 
is consistent with current LHC data [13j. 

In our analysis, we consider the LHC experiment with the center of mass energy 
of ^/s = 14TeV. Then, the signal cross section, which is the sum of the cross sections 
for gluino and squark productions, is estimated to be 

a~gg = 0.129pb, aqg = 0.922pb, aqq = 0.879pb. (7) 

We focus on the decay chain involving a right-handed squark, as shown in Fig. [H As 
a result, a typical signal event consists of two energetic jets, four r-leptons (r-jets or 
leptons), and a missing energy in the transverse direction. 

For parton-level event generation and hadronization, we employ the HERWIG 
code [T8l[T9]. Generated events are passed through the PGS code [20] for simulat- 
ing detector effects. Fake r-jets from QCD processes and heavy meson decays are 
involved in the study. For tau-jets, we smear the transverse vertex position of the 
parton using Gaussian distribution with the error Adj = 10/im. 
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Mass (GeV) Branching fractions 



1096.6 Br(^ qq) = 0.89. 

951.1 Bt{ul xfd) = 0.34, Bt{ul Xid) = 0.32, 
Br{uL xlu) = 0.18, Br{uL x'^u) = 0.15. 

922.0 Bt{ur x^u) = 0.96. 

197.3 Br(x? ^ f^r"^) = 0.35, Br(x? ^ e^e^) = 0.32, 

Br(x? ^ = 0.32. 

130.0 Br(ei? f^r'^e) ~ 1.00. 

126.2 Br(fi ^ tG) = 1.00. 



Table 1: Masses and branching fractions of sparticle in our representative point. 



3.2 Charge subtraction method 

Background reduction is the most important task in our analysis to determine the 
mass spectrum of sparticles, because all kinematical endpoints do not have sharp 
edge structures due to the energy leakage by Ur emissions from r-decays. Expected 
backgrounds in those measurements are as follows: (i) A number of fake r-jets are 
expected at the hadron collider, (ii) A signal event results in multiple r-leptons, and 
hence there are combinatorial backgrounds in the analysis involving r-jets. 

In order to reduce these backgrounds, we adopt the method of the charge subtrac- 
tion. We expect four r-leptons (two r+ and two r^) in one event. We therefore have 
three ways to pair the r leptons, (r^, r^)i, (r^, r^)2, and (r^, r^). In each event (in- 
volving four rs), we take the data using the method {r^, r^)i -|- {r^, r^)2 — {t^, r^), 
then the wrong opposite-sign pairing is expected be canceled by the subtraction of 
the same-sign pairing. This method works very well when r-leptons are produced in 
the process x° — >■ f — >■ G; this is due to the fact that x° decays into f+r~ and f~T^ 
with equal probability. 

In addition, this method can be applied to the determination of the M^yCnear)- 
endpoint. We simply collect {j, r, r) events using the charge subtraction. In each 
paring, r'^^^^'^^ is identified as the r lepton which has a shorter impact parameter. It 
should be also noted that the charge subtraction method can reduce backgrounds 
from fake r-jets from QCD processes, because the QCD fake events are charge-bhnd 
in a good approximation at high energy processes. 
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Selection cut 


SUSY 


tt 


(0) Generated events 


273,600 


49,610,000 


(1) i)T > ISOGeV 


190,504 


1,644,160 


(2) # of leptons > 4 


70,641 


55,537 


(3) # of b-jets = 


56,228 


41,673 


(4) ji with pt > 200GeV 


49,819 


5,915 


(5) j2 with Pt > lOOGeV 


41,007 


1,984 



Table 2: Numbers of signal (SUSY) and background (tt) events after applying kinematical 
cuts with C = lOOfb^^. All SUSY processes are included in our event generation. 

3.3 Kinematical cuts to reduce ti backgrounds 

The most serious SM background for our study is the ti production. Thus, we 
concentrate on this background. In order to reduce this background, we impose 
following kinematical cuts: 

• Large missing transverse energy, > 150GeV. 

• At least, four leptons, e, /i with pt > 20GeV or r-jet with pt > 25GeV. 

• Two hard jets, ji(j2) with px > 200(100)GeV (and no b-jets). 

Here, we use the label, i = 1 or 2, for a jet (ji) in decreasing order of pt- We also take 
account of both electrons and muons, because r-lepton often decays leptonically and 
selectron and smuon decay into gravitino by emitting electron and muon directly. 
The requirement for jets to have pt > 200(100)GeV is very important to reduce the 
tt background, because jets with px > mw {m^y is the mass of weak gauge boson) 
from top quark decays are rather rare. In Tabled we summarize the cut flow in our 
simulation study with assuming that the integrated luminosity is lOOfb^^, where we 
take the gravitino mass of 9.7eV again. 

It is possible to apply tighter kinematical cuts for further reductions of the back- 
grounds. For instance, the requirement > 200GeV, ptOi) > 250GeV and pt(j2) > 
150GeV in addition to the basic kinematical cuts shown in Table [2] will reduce 80% 
of the ti background (we have, as a result, ~400 ti events), while this also reduces 
30% of SUSY signals (we have, as a result, ~28000 SUSY events). 
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4 Simulation results 



We are now in position to present several results of our simulation study, which 
are obtained based on arguments in previous sections. We first show the results for 
measurements of sparticle masses, and discuss how accurately these masses can be 
determined at the LHC We next show that the lifetime of the stau may be deter- 
mined by using the distribution of the transverse impact parameter of r-jets from 
the stau NLSP decay. We estimate how accurately the lifetime can be determined. 

4.1 Sparticle masses 

The strategy to determine the mass spectrum is the use of kinematical endpoints 
of several variables. We study how the endpoints behaves using generated events 
which pass through the basic cuts discussed in section 13.31 For the simulation study 
of sparticle mass measurement, the decay length of the stau NLSP is set to be 500/im 
(corresponding to the gravitino mass of 9.7eV). 

4.1.1 Endpoint on Mrr 

The first kinematical variable used in the analysis for the mass spectrum is the in- 
variant mass of two r-leptons in the decay chain of a squark. In Fig. |2] (upper panel), 
the distribution of the invariant mass (M^^) after applying the charge subtraction 
method is shown. We can see a clear edge at Mrr — 150 GeV. In order to extract 
the location of the endpoint, we use the following fitting function. 



where A, B, M^*, and C are parameters to fit the shape of the distribution around 
the endpoint. With the use of this bilinear function for the fitting, the location of 
the endpoint is determined to be M^* = 151.2 ± 14. 5GeV. Notice that the underlying 
value (the input value on the simulation) is 151.6 GeV. 

4.1.2 Endpoint on M^^ I near) 

Second kinematical variable we use is the invariant mass between r-lepton and jet 
emitted by the decay of a squark. Using information about the impact parameter, it 




11 




J?/ndf 
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B 
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167.07532±100. 57011 
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Figure 2: (Upper panel) Distribution of the invariant mass between two tau-jets, Mrr- 
(Middle panel) Distribution of the invariant mass between hard jet and near r-jet, Mj^(nca.v) ■ 
(Lower panel) Distribution of the MT2,jj variable defined in Eq. ([6]) with m^iss being zero. 
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is possible to distinguish near and far r-leptons with high efficiency. For each (r, r)- 
pair, we identify the r-jet whose track has a larger impact parameter than the other 
as the far tau-jet, while the r-jet with a smaller impact parameter is regarded as the 
near r-jet. After this identification, both combinations of (ji, t^^^^^")) and {j2, t^^'^^^^) 
are used to calculate Mq^(noar). In the analysis, we also require that the pair of two 
tau-jets should satisfy M^-t < M^* = 155.72GeV in order to reduce fake-QCD and 
combinatorial backgrounds. 

The distribution of the invariant mass (Mg^(near)) after applying the charge sub- 
traction method is shown in Fig. [2] (middle panel). The endpoint is, again, fitted by 
using the bilinear function given in Eq. ([8]). It then turns out that the location of 
the endpoint on Mg^( near) IS M*^^\„,,,) = 700.0 ± O.lGeV. Notice that the underlying 
value is now 692.3GeV. 



4.1.3 Endpoint on Mt2. 



n 



The last kinematical variable is MT2,jj defined by two hard jets as visible particles, 
as mentioned in section [2^21 In our analysis, we take the test mass (mmiss) in Eq. 
to be zero, so that the endpoint of this kinematical variable gives M^2^j{0) = (m? — 
m'io)/mq. The distribution of Mr2,ji(0) is shown in Fig. [2] (lower panel). It can 
be seen that a very clear endpoint exists at Mr2jj(0) — 870GeV. As in the cases 
of previous kinematical variables, we fit the shape of the distribution around the 
endpoint by the bilinear function. The endpoint of the distribution is then obtained 
as M^\ -.{<d) = 875.1 ± 6.7GeV (the underlying value is 879.8GeV). 

4.1.4 Mass determination 



When the gravitino mass is neglected, the masses of sparticles, rriq, m^o, and mf^, 
are determined by three kinematical endpoints of the variables, M^^^, M™('^g^,), and 
M™2^^j(0). Using analytic expressions for the these endpoints shown in eqs.(l3]), (jl]), 
and ([2]), the masses of sparticles are determined by minimizing the following 
function, 



Am 



]Vf max _ j^&t -1 2 

TT TT 



+ 



/If max _ fit 

qT'- 



-\ 2 



+ 



Mf-(0)-M|* ,(0) 



rfit 

AMfi*_^.^.(0) 



2 



(9) 



^(near) 

where M^* denotes the center value of the measured endpoint, and AM*^* is its 
(statistical) error. After minimizing by varying the input values, m^, rn^o^ and 
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Figure 3: Distribution of the transverse impact parameter of r-jets after applying the 
basic kinematical cuts in previous section. We show distributions in the SUSY model with 
CTf = 1, 100, 500 and 900/im, respectively, and also show the distribution in tt production. 

mf^, we obtain following results; the right-handed squark mass is rriq = 915.9 ± 
6.4GeV (the true value is 922.0GeV), the lightest neutralino mass is m^o = 193.4 ± 
19.5GeV (the true value is 197.3GeV), and the lightest stau mass is m^^ = 120.5 ± 
IS.lGeV (the true value is 126.2GeV). Here, the gravitino mass is taken to be zero 
(i.e., negligibly small). 



The strategy to determine the lifetime of the stau NLSP is the use of the distribution 
of the transverse impact parameter. After showing the distribution for several input 
values of CTf, we discuss how accurately the lifetime can be determined at the LHC. 

4.2.1 Distribution of the impact parameter 

Distribution of the transverse impact parameter (dj) obtained from hadronically 
decays of r-leptons is shown in Fig. [3] with the use of generated events which are 
passed through the kinematical cuts discussed in previous section. Four distributions 
are shown in this figure with choices of the decay length of the stau NLSP to be CTf = 
1, 100, 500, and 900/im, respectively. It is clearly seen that the distribution depends 
on the decay length of the stau NLSP as expected. 



4.2 Lifetime of the stau NLSP 
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It is also seen that a number of events are found in inner bins with dj ^ CTf and 
that a broad tail-structure exists in the region dj ^ cr^ . In fact, in both regions, 
backgrounds are expected to contribute to the distribution. In small dj region, the 
distribution is dominated by background r-jets such as Q CD-originated fake ones. 
Though most of those backgrounds do not have finite dj at the parton-level, the 
backgrounds acquire finite values of dj at the detector-level because of the limited 
resolution for the vertexing. On the other hand, in the region of dj ^ cr^ , back- 
grounds come from decays of heavy hadrons such as D ot B mesons. Those hadrons 
sometimes produce fake r-jets after flying a sizable distance. In order to eliminated 
these backgrounds efficiently, we vary and optimize the upper and lower endpoints 
of the bins which are used for the analysis as we change cTf"^^^^ (where ct^^'^^^'' is 
the test value of the decay length used to generate a template of dj distribution). 

4.2.2 Lifetime estimation 

According to the strategy to estimate the lifetime discussed in section 12. 3[ we now 
study how well we can constrain the lifetime of fi. First, we use the template 
generated with the true value of the averaged velocity of the stau NLSP, which is 

= 0.88c in our representative point. In order to see how the result depends on the 
underlying value of the lifetime of fi, here we use several values of CTf in generating 
the impact-parameter distribution. In addition, in preparing the templates for dj 
distributions, TAUOLA library [21] is used to simulate the r decay event, which 
enable us to deal with chirality and finite lifetime of r-leptons. The range of the test 
value is taken to be lOyum-llOO/im every lO/xm. 

With the use of the templates, we perform x^-analysis to determine the lifetime. 
In our analysis, only r-jets satisfying 0.5 x ct-^'^'^^^ < dj < 2.0 x cr^*'^'^*^ with small 
pseudo-rapidity, \ri\ < 1.0, are adopted. The size of bin used in each distribution is 
set to be 0.15 x cr^*'^'^*^ and, as a result, we have 10 bins in total. The variable to 
estimate the lifetime of the stau NLSP is therefore given by 



where A^^ (cr^ ) denotes the number of signals in the i-th bin obtained by the 
template for a given ct^^'^'^^\ while N^'^^^'^ is the one obtained by using generated 




(10) 
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Figure 4: x^'Values as a function of tlie test lifetime (the test decay length) of the stau 
NLSP, CT^°^^\ for the underlying values of CTf = 100 (red, solid), 300 (green, dashed), 500 
(blue, dotted), 700 (violet, dot-dashed), and 900 (cyan, dot-dot-dashed) /xm. Here, the 
averaged velocity is taken to be i3f = 0.88. 



events. We only involves the statistical error as ANi = y N^^^'^ in the analysis. The 
degrees of the freedom in this x^-test is therefore (10 — 1) =9, and the hypothesis 
is excluded at 95% C.L. when > 16.92. 

Resultant x^-values as a function of the test lifetime (the test decay length) of 
the stau NLSP, cTf°^^\ for underlying values of CTf = 100 (red, solid), 300 (green, 
dashed), 500 (blue, dotted), 700 (violet, dot-dashed), and 900 (cyan, dot-dot-dashed) 
/im are shown in Fig. |H From these results, the lifetime (the decay length) of the 
stau NLSP in each case is determined at 95% C.L. to be 



We can see that, if the correct value of /3f is used, the analysis based on the impact 
parameter distribution gives a good estimate of the lifetime (the decay length) of 
the stau NLSP with accuracy of about 30% when cTf > 100/im. If more precise 
information about the velocity distribution of fi is available, better estimate of the 
lifetime may be obtained. 




50/im < CTf < llO/xm (underlying CTf 

240/im < CTf < 330/im (underlying CTf 

410/xm < CTf < 540/im (underlying CTf 

570/im < CTf < 800/im (underlying CTf 

810/xm < CTf < 1060/im (underlying CTf 



100/im), 
300/-fm), 
500/im), 
700/im), 
900/-im). 



(11) 
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So far, we have neglected the uncertainty arising from the determination of the 
velocity distribution of fi. As we have mentioned, the detailed study of the uncer- 
tainty in the velocity distribution is beyond the scope of this paper. In our analysis, 
however, we estimated the uncertainty of the Pf determination related to the errors 
in mass measurements and also to the production process in order to demonstrate 
that the Pf can be obtained with some accuracy. 

First, in order to study the effects of the errors in the mass measurements, we 
generated the events using the different sparticle mass spectrum from our represen- 
tative point. Here, we used the mass spectrum predicted from the simple gauge 
mediation model except for fi because we found that the error of is the largest 
among the reconstructed masses in the previous subsection. Then we generated the 
full SUSY events by varying m^^ by ±20GeV, and found that the value of (3f changes 
by ~ 0.01. In addition, the dominant SUSY process may not be well understood 
in the actual situation. If so, it may be reasonable to estimate /3f by assuming the 
process we use for our analysis, which is the process shown in Fig. [H We generated 
events corresponding to such a process (using the correct mass relation). Then, is 
found to be ~ 0.93. Thus, a relatively larger uncertainty of A/3f ~ 0.05 is expected 
if the dominant SUSY process cannot be understood. 

To see how this affects the determination of the lifetime, we calculate the 
variable using the template with jSf = 0.83 and 0.93. The results are shown in 
Figs. [5] and [61 We can see that error related to the uncertainty to the averaged 
velocity is ~ 50 — lOO/xm if APf ~ 0.05. Even with such an uncertainty, we can still 
have a relatively good determination of Tf. Thus, the impact parameter will be a 
powerful tool to measure the lifetime at the LHC 

Finally, we consider how well we can estimate the gravitino mass if we assume 
that fi decays into the gravitino and r-lepton. In the case with the underlying 
gravitino mass of 9.7eV, we have shown that the lifetime of the stau is estimated 
between 410/im and 540/im for f3f = 0.88. In addition, the uncertainty related to f3f 
is estimated to be about 40/im, as one can see from Figs. [5] and [61 Then, by using 
the measured stau mass rrif = 120.5 ± 18. IGeV and the center value of the estimated 
NLSP lifetime cr^ = 475/im, we obtain m3/2 = 8.3 ± 1.8 ± 0.8 ± l.OeV, where the 
errors originate in the uncertainties of m^, cTf and (3^, respectively. 
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Figure 5: Same as Fig. HI except for f3f = 0.83. 
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Figure 6: Same as Fig. [U except for f3f = 0.93. 



5 Summary 

We have proposed a method to determine the mass spectrum of sparticles and the 
hfetime of the stau NLSP at the LHC when the mass of the gravitino LSP is of the 
order of lOeV. Though the decay length of the stau NLSP is very short, which is of 
the order of 100-1000/im, it is still possible to deeply study the model by utilizing 
the transverse impact parameter of tau-jets from the decay of the stau NLSP. 

We have first discussed the mass measurement of sparticles using a typical cas- 
cade decay chain of a squark shown in Fig. [H This SUSY event involves, at least, 
four r-leptons, which makes it difficult to analyze the signal because of combinato- 
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rial backgrounds. Information about the impact parameter of the tau-jet, however, 
resolves the problem, and we have shown that the mass spectrum of sparticles can 
be determined accurately through the kinematical endpoints of Mtt, Mj^(nea.T), and 

MT2,jj. 

We have also discussed the determination of the lifetime of the stau NLSP using 
the distribution of the transverse impact parameter. The impact parameter depends 
not only on the lifetime but also on the velocity of the stau NLSP. We have therefore 
developed a strategy to estimate the velocity by utilizing a simulation with infor- 
mation about the mass spectrum obtained in the previous stage. We have shown 
that, if the velocity distribution of fi is somehow understood, the lifetime of the 
stau NLSP is determined with the accuracy of about 30% as far as its decay length 
is larger than ~ 100//m. Thus, if the underlying model of the SUSY breaking is 
low-scale gauge mediation with the gravitino mass of O(10)eV, the LHC may have 
a chance to acquire some information about the gravitino mass. 
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